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Abstract. It has been suggested that relativistic blast waves 
may power the jets of AGN and gamma-ray bursts (GRB). 
We address the important issue how the kinetic energy of col- 
limated blast waves is converted into radiation. It is shown 
that swept-up ambient matter is quickly isotropised in the blast 
wave frame by a relativistic two-stream instability, which pro- 
vides relativistic particles in the jet without invoking any ac- 
celeration process. The fate of the blast wave and the spectral 
evolution of the emission of the energetic particles is therefore 
solely determined by the initial conditions. We compare our 
model with existing multiwavelength data of AGN and find re- 
markable agreement. 
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1. Introduction 

Shortly after the detection of more than 60 blazar-type active 
galaxies as emitters of MeV-GeV 7-rays by EGRET on board 
the Compto n Gam ma-Ray Observa tory ( CGRO) (von Mon- 
tigny et al. 



1995 



1997), observations with 



Mukherjee et al. 
ground-based Cerenkov telescopes have shown that at least for 
BL Lacertae objects the 7-ray spectrum can be tr aced t o more 
than a TeV observed phot on en ergy (Punch et al. 1992 Quinn 
et al. |1996|, Catanese et al. |l998j Rrennrich et al. |1997|). In these 



sources the bulk of luminosity is often emitted in the form of 
7-rays. The typical 7-ray spectra of blazars are well described 
by power laws in the energy range between a few MeV and a 
few GeV (Mukherjee et al. 1997), but display a break around an 
MeV (McNaron-Brown et al. 1995 ). At GeV energies the spec- 
tra show a tendency of getting harder during outbursts (Pohl et 
al. 1997). The 7-ray emission of the typical blazar is highly 



variable on all timescales down to the observational limits of 
days at GeV energies and hours at TeV energies (e.g. Mattox 
et al. 1997, Gaidos et al. 1996 1. Large amplitude variability 



at TeV energies is generally not associated with correspond- 



ing variability at GeV energies (Lin et al. 1994, Buckley et al 



1996) 



Correlated monitoring of the sources at lower frequencies 
has demonstrated that 7-ray outbursts are accompanied by ac- 
tivity in the optical and radio band (Reich et al. 1993, Miicke et 
al. 1996, Wagner 1996[ ). Emission constraints such as the com- 
pactness limit and the Elliot-Shapiro relation (Elliot & Shapiro 
|l974t ) are violated in some 7-ray blazars, even when these lim 



its are calculated in the Klein-Nishina limit (Pohl et al. 1995 



Dermer & Gehrels 1995), which implies relativistic bulk mo- 



tion within the sources. This conclusion is further supported 
by observations of apparent superluminal motion in many 7- 



Kingham 1997a, 1997b) 



ray blazars (e.g. Pohl et al. 1995 Barthel et al. 1995, Piner & 



Published models for the 7-ray emission of blazars are usu- 
ally based on inverse Compton scattering of soft target pho- 
tons by highly relativistic electrons in the jets of these sources. 
The target photons may come directly from an accretion disk 



(Dermer et al. 1992 Dermer & Schlickeiser 1993), or may 



be rescattered accretions disk emission (Sikora et al. 1994), 



or may be produced in the jet itself via synchrotron radiation 



(e.g. Bloom & Marscher 1996 and references therein). In all 



the above models the jet and its environment is assumed to 
be optically thin for 7-rays in the MeV-GeV range. A com- 
mon property of the inverse Compton models is the empha- 
sis on the radiation process and the temporal evolution of the 
electron spectrum and the neglect of the problem of electron 
injection and a cceleration (for an exception see Blandford & 
Levinson 1995 ). The modeling of the multifrequency spectra 
of blazars often requires a low energy cut-off in the electron in- 
jection spectrum, but no high energy cut-off. The usual shock 
or stochastic electron acceleration processes would have to be 
very fast to compete efficiently with the radiative losses at high 
electron energies, and a high energy cut-off should occur under 
realistic conditions, but no low energy cut-off (e.g. Schlickeiser 



1984), in contrast to the requirements of the spectral modeling. 



When efficient, but possibly slower, acceleration of protons 
is assumed, photomeson production on ambient target photons 
can provide many secondary electrons and positrons with a 
spectrum a s requ ired by the mu ltifreq uency modeling (Kazanas 
& Ellison I1986L Sikora et al. I1987L Mannheim & Biermann 



1992). Such systems are usually optically thick and a pair cas- 
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cade develops. In these models multi-TeV 7-ray emission can 
be easily produced, but the observed short variability timescale 
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places extreme constraints on the magnetic field strength, be- 
cause the proton gyroradius has to be much smaller than the 
system itself, and on the Doppler factor, because the intrinsic 
timescale for switching off the cascade is linked to the observed 
soft photon flux via the energy loss rate for photomeson pro- 
duction. For the most rapid 7-ray outburst of Mrk 421 (Gaidos 
et al. 119961) this leads to B > 10 G and D > 100. 



Another class of models features hadronic collisions of a 
collimated proton jet with BLR cl ouds e ntering the jet (e.g. 
Dar & Laor ll997| ; Beall & Bednarek[l99§). These models have 
two theoretical difficulties: the proton beam is weak compared 
with the background plasma and therefore quickly stopped by 
a two-stream instability. Also the BLR clouds are usually op- 
tically thick and thus the efficiency of the system is drastically 
reduced. 

In this paper we consider a strong electron-proton beam that 
sweeps up ambient matter and thus becomes energized. The ba- 
sic scenario is similar to the blast wave model for 7-ray bursts 
(GRBs) which successfuly explains the time dependence of the 
X-ray, optical and radio afterglows (e.g. Wijers et al. 1997, 
Vietri 1997a, Waxman 1997a). There the apparent release of 
~ 10 52 ~^2~erg of energy in a small volume leads to the forma- 
tion of a relativistically expanding pair fireball that transforms 
most of the explosion energy into kinetic energy of baryons in 
a relativistic blast wave. We assume that a similar generation 
process also powers the relativistic outflows in active galac- 
tic nuclei (AGN), but that the outflow is not spherically sym- 
metric and highly channelled along magnetic flux tubes into a 
small fraction of the full solid angle due to the structure of the 
medium surrounding the point of energy release. 

We shall address the important issue how the kinetic energy 
of such channelled relativistic blast waves is converted into ra- 
diation. Existing radiation modelling of GRBs and AGNs (see 
1997bt Waxman |1997b| ; Bottcher & Dermer |1998| ) 



e-g 



Vietri 



are very unspecific on this crucial point. Typically it is as 
sumed that a fraction £ of the energy in nonthermal baryons 
in the blast wave region at position x is transformed into a 
power-law distribution of ultrarelativistic cosmic ray protons 
with Lorentz factors T(x) < 7cr < Ymax m the fluid frame 
comoving with a small element of the blast wave region that 
travels with the bulk Lorentz factor T(x) — T (x/xo)~ 9 after 
the deceleration radius xq = 2.6 • 10 16 (E52 / noT^ 30 o) 1 ^ 3 cm 
through the surrounding mediu m of density no. For elect rons i t 
is argued (Katz & Piran 1997 , Panaitescu & Meszaros 1998 ) 
that their minimum Lorentz factor is 7 c . m i n = (m p /m e )T 
since they are in energy equilibrium with the protons. Here 
we investigate this transfer mechanism from the channelled 
blast wave to relativistic protons and electrons in more de- 
tail. In Sect. 2 we consider the penetration of a blast wave 
consisting of cold protons and electrons with density rib with 
the surrounding "interstellar" medium consisting also of cold 
protons and electrons on the basis of a two-stream instability. 
Viewed from the coordinate system comoving with the blast 
wave, the interstellar protons and electrons represent a proton- 
electron beam propagating with the relativistic speed V(x) = 
— c(l — r _2 (.x)) 1 / 2 antiparallel to the x-axis. We examine 



the stability of this beam assuming that the background mag- 
netic field is uniform and directed along the x-axis. We demon- 
strate that very quickly the beam excites low-frequency magne- 
tohydrodynamic plasma waves, mainly Alfven-ion-cyclotron 
and Alfven- Whistler waves. These plasma waves quasi-linearly 
isotropise the incoming interstellar protons and electrons in the 
blast wave plasma. In Sect. 3 we investigate the interaction pro- 
cesses of these isotropised protons and electrons (hereafter re- 
ferred to as primary protons and electrons), which have the rel- 
ativistic Lorentz factor T(x), with the blast wave protons and 
electrons. Since the primary protons carry much more momen- 
tum than the primary electrons, inelastic collisions between 
primary protons and the blast wave protons generate neutral 
and charged pions which decay into gamma rays, secondary 
electrons, positrons and neutrinos. Both, the radiation prod- 
ucts from these interactions, and the resulting cooling of the 
primary particles in the blast wave plasma, are calculated. By 
transforming to the laboratory frame we calculate the time evo- 
lution of the emitted multiwavelength spectrum for an outside 
observer under different viewing angles. Momentum conser- 
vation leads to a deceleration of the blast wave that is taken 
into account self-consistently. Since we do not consider any 
re-acceleration of particles in the blast wave, the evolution of 
particles and the blast wave is completely determined by the 
initial conditions. 

2. Two-stream instability of a proton-electron beam 

2.1. Physical model and basic equations 

As sketched in Fig. 1 we consider in the laboratory frame (all 
physical quantities in this system are indexed with *) the cold 
blast wave electron-proton plasma of density n£ and thickness 
d* in x-direction running into the cold interstellar medium of 
density ng, consisting also of electrons and protons, parallel 
to the uniform magnetic field of strength B. In the comoving 
frame the total phase space distribution function of the plasma 
in the blast wave region at the start thus is 

/>*.* = 0)= in^( F l)j(pi) 



2np* ± 
1 

2irp*. 



nlS(p* ± )6(pf, - P) 



(1) 



where 



P = YmV = Tmfic 



rac 



v 7 ^ 



1) 



The Lorentz transformations (e.g. Hagedorn 1963 , p.41) to the 
blast wave rest frame with momentum and energy variables 

(P±,P\\,E) are 



pl =p±, pI = r(pii 



c 



E* =T(E+ (3cp {] ) . (2) 



Using the invariance of the phase space densities we obtain 
for the phase space density in the blast wave rest frame 



f(p,t = 0) = f*(p*(p),t = 0) 

ni6(p±)S(p\\ - P) n h 8(pi_)8{p\\ 
2irp± 2irp± 



(3) 
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Fig. 1. Sketch of the basic geometry. The thickness of the chan- 
neled blast wave d, measured in its rest frame, is much smaller 
than its halfdiameter. The blast wave moves with a bulk Lorentz 
factor T through ambient matter of density n\. 



where the number densities are 

ni^TriQ, n h = n* h /Y. 

The blast wave density transformation is inverse to that for the 
density of interstellar matter because of the different rest sys- 
tems. 

In terms of spherical momentum coordinates pu — p[i and 
p± = P\J\- — y, 2 the distribution function (Q) reads 

1 „ / W , „s 1 



f(p,t = 0) = 



niS(fi + l)S(p - P 



27rp 2 

n bfh(p, t 



Anp 2 
0) + riifi(p,iJ,,t) , 



n h S (p) 



(4) 



where /i = pii /p is the cosine of the pitch-angle of the par- 
ticles in the magnetic field B pointing to positive z-direction. 
Obviously, the incoming interstellar protons and electrons are a 
beam propagating antiparallel to the magnetic field direction in 
the blast wave plasma. For ease of exposition we assume that 
the blast wave plasma is cold, i.e. /b oc S(p). 

The blast wave density n b = 10 8 n£ g is much larger than 
the interstellar gas density = In*. We want to calculate the 
time tpi it takes to isotropise the incoming interstellar protons 
and electrons; if this relaxation time is much smaller than d/c 
an isotropic distribution of primary protons and electrons is ef- 
fectively generated. 

If the beam distribution function fa in Eq. ^ is unstable in 
the blast wave plasma it will excite plasma waves. Therefore 



we have to find the dispersion relation u> 
perturbations of the form 

5E = SE exp(i(fc • x — u;t)) . 



w(k) which are 



(5) 



It has been shown (Kennel & Wong 1967 Tademaru 1969 ) that 
waves propagating obliquely to the magnetic field are thermally 
damped by the background (=blast wave) plasma, so that we 
concentrate on parallel (fc||.B) waves. One finds two implicit 
dispersion relations (Baldwin et al. 1969; Achatz et al. 1990), 
one for electrostatic waves (SE\ \B) 



n 2 2 2 

= LU - UJ p p - CJ p o 



1 



n h T 3 

and one for electromagnetic (5E _L B) 



U! 2 LU 2 
p,c 



(6) 



J+kV] 



j-O c 



= 



(7) 



PJ — y Airribqf /mj denotes the plasma frequency and 



where to 

flj = qjB/rrijC the nonrelativistic gyrofrequency of particle 
species j. In our notation negative (positive ) frequencies wr. 
refer to right-handed (left-handed) circularly polarised waves, 
respectively. Moreover, positive phase speeds (wR,/fc > 0) in- 
dicate forward propagating waves and negative phase speeds 
(ojYi/k < 0) backward propagating waves. Here we will con- 
centrate on the excitation of the electromagnetic waves. The 
excitation of electrostatic turbulence and its backreaction on 
the particle distribution function will be discussed in a forth- 
coming publication. Because 



rii/nb = T 2 n* /n h = 10 4 ri 00 n*/n b .8 << 1 



(8) 



the beam is weak and the contributions from the beam in Eq.(Q) 
can be described as perturbations to the dispersion relation in 
a single component electron-proton plasma (?i; = 0). Its solu- 
tions 



(WR 



n P )a 



are 

(a) Alfven waves 



2;„2 



(9) 



(10) 



at frequencies |wr| << 51 p , where Va = -B(47rm p nb)" 1 / 2 is 
the Alfven speed and cjr = dt(ui). The dispersion relation ( |Io| ) 
accounts for four types of Alfven waves: forward and backward 
moving, right-handed and left-handed polarised; 
(b) Whistler waves 



ujr, = Sl e k 2 c 2 /u) 



p.c 



(11) 



at frequencies between Q c < < — Q p . The dispersion 
relation (|llj) describes right-handed polarised waves (because 
il c < 0) that propagate forward for negative k < and back- 
ward for positive k > 0. 
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All of these are stable (tp = ?y(u>) = 0) if the beam par- 
ticles are absent. As is explained in Achatz et al. (1990 1 the 
beam protons and electrons, which under the given weak-beam 
condition do not affect each other, are however able to trigger 
each its instability which occurs if the resonance conditions 



uj R = -kV + (Op/T), uj k = -W+ (fia/r), 



(12) 



can be satisfied. 

The time-dependent behaviour of the in tensities J(fc, t) of 
the excited waves is given by (Lerche 1967, Lee & Ip 1987) 



dl u 

where the growth rate ip is 



(13) 



oo E 2 
dE 



8uj k 
1 - 

1 ^ N 



sgn(k)J]a;2 i(mic) 3 

i 

2 df t 



dfi 



5{fx 



), (14) 



where N = ckj lj r is the index of refraction and further Ei = 
\/l + xf whith Xi = Qifi/kc. 

To describe the influence of these excited waves on the 
beam particles we use the quasilinear Fokker-Planck equation 
(e.g. Schlickeiser 1989) for the resonant wave-particle interac- 
tion. For Alfven waves and for Whistler waves the index of 
refraction 



\Na\ 



1 



( c /y A ) 2 » i 



in 



w 



(15) 
(16) 



is large compared to unity, so that the Lorentz force associated 
with the magnetic field of the waves is much larger than the 
force associated with the electric field, so that on the shortest 
time scale these waves scatter the particles in pitch angle \i but 
conserve their energy, i.e. they isotropise the beam particles. 
The Fokker-Planck equation for the phase space density then 
reads 



df, 



dt dfj} 



dfi- 

"» du- 



al) 



where the pitch angle Fokker-Planck coefficient is determined 
by the wave intensities I n 



D„ 



V 7rOf(l-M 2 ) 
^ 2B 2 



dk I n (k) 6(u> R — kvfi — Qi) 



Using the dispersion relation of Alfven waves (uj r ~ Vkk), Eq. 
(18) indicates that beam protons and electrons resonate with 
waves at wavenumbers given by the inverse of their Larmor 
radii times \x, k ~ Cli/Fvfi — (/iRi)^ 1 . For protons these 
wavenumbers correspond to Alfven waves. For electrons these 
wavenumbers correspond to Alfven waves if the bulk Lorentz 



factor is above | fi\(m p /m c ) = 1836|/i|, and Whistler waves for 
smaller Lorentz factors. We concentrate here on the isotropisa- 
tion by Alfven waves mainly for two reasons: 

(1) the bulk of the momentum of the inflowing interstellar par- 
ticles is carried by the protons so that they are energetically 
more important than the electrons; 

(2) for Lorentz factors V > |/i|m p /m e the isotropisation of 
electrons is also caused by scattering with Alfven waves; for 
smaller Lorentzfactors the mistake one makes in representing 
the Whistler dispersion relation still by the Alfven dispersion 
relation is relatively small. 

2.2. Self excited Alfven waves 
For Alfven waves (|l0|) we obtain 

dJ(iv R ) 



dujv 



2lj r c*/V£ = 2c 2 k/V A 



(19) 



which is positive for forward (k > 0) moving waves and neg- 
ative for backward (k < 0) moving waves. According to Eq. 
( p"4| ) we obtain for the growth rate of forward (+) and backward 
(-) moving Alfven waves 



with 



(20) 



(21) 



Now it is convenient to introduce the normalised phase space 
distribution function of the beam particles 



/i(p,M»*) 



s(e - r) 



27r(m iC ) 3 r(r 2 - 1)1/3 



(22) 



where E = yl + (p/mc) 2 , so that Eq. (gTj) reduces to 

f^E^^r-v/TT^ni-A] x 



^(M-**(r 2 -i)-V3) 



(23) 



where H[x] = 1(0) for x > (<)0 is the step function. Sum- 
ming over protons and electrons separately, and introducing 
the proton and electron Larmor radii, R p = c\/T 2 — 1/0 P 



R e = cy/T 2 — l/|f2 c |, respectively, Eq. (23 ) becomes 

^ = K|fc|)Ex=e,p 



£[1 - (kR x )- 2 ]H[\k\ - R- 1 }^^ + (Rxk)- 1 ] 



(18) with 



6(1*1) 



vr Va rii uj v 



2 c n\yT c\k\ 
According to Eqs. ( |l3| ) and (p 

01+ (t) 



• written in the form 



at 



(24) 



(25) 



(26) 
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^- = -^I-(t) (27) and 



we derive 



£ / dtI ± (\k\,t)-^5( f i-(R c k)- 1 ) = 



d(I+(t)I-(t)) _ n (2g) , 1/Refe 



± 

77/ ~" a-WEck]- 1 ) I dn[F e (n,t) - F c (fi,t = 0)] (38) 



and 

d[I + (t)-I.(t)\ 



Eqs. (I37p and (B8h can be inserted in Eq. (B2h yielding 



„ np rprr fl ~ (fc-Rp) ] 



m p a p (\R p k\ 1 ) 

-sgn(f2 x n) 



V H[\k\ - R- 1 } I Rx " dft[F x (n,t)-F x (ii,0)] (39) 



or when integrated, 

J+ (t) J_ (t) = /+(* = 0)1- (t = 0) (30) 
and 

[/+(«) - /_(«)] - [7+(i = 0) - /_(* = 0)] = Z(Jfe) (31) The system of the coupled equations ©flj|p) and © de- 
scribes the temporal development of the particle distribution 

with functions F p und F c under the influence of self-excited Alfven 

Z {k) = b(\k\) Y ^[1 - (kR^}H[\k\ i?- 1 ] x ™ v< ® Propagatingeitherforward( + )orbackward(-).Her e the 

m ^ boundary conditions are as follows: in the beginning it = 0) 

t x c p there is the mono-energetic beam distribution (jlj), i.e. in terms 

j dt ^S{fi + (R^k)- 1 ) (l+(t) + I- (t )) (32) of the normalised distribution @ 

F Pie (jM,t = 0) = 6(ji+l) (40) 

which relates the wave intensities of forward (+) and backward 

(-) moving waves at any time t to the initial wave intensities at The final state of the isotropisation phase is reached at time Ti 

time t = 0. when both growth rate and temporal derivative of the distribu- 

In terms of the normalised distribution functions (|||) the tion disappear, i.e. when dFi/dfi = 0. Consequently 
proton and electron distribution functions evolve according to 

Eqs.©-©, " F P)e (/i,t = T0 = i (41) 

— - — — [a p (|/x|) I±(|i?p/x| _1 ) — -] (33) At this time the Alfven waves have completely isotropised the 

^ ^ ± ^ beam distrution. In order to estimate this time scale Tj or the 

associated isotropisation length A we consider the extreme case 

F)F f) 3F 

— _ [ ac (|^|) \^ Jj- (| Jge/xp 1 ) j (34) that the wave spectrum is constant and equals the wave inten- 
ds dfi ± dfj, sity spectrum after the isotropisation. By this method an ap- 
proximation of the pitch angle Fokker-Planck coefficient D MM 
with is possible so that one gets a strict lower limit to the isotropisa- 
tt v (i _ ^\ it v(l — /j 2 ) ^ on l en S m A. By demonstrating that this length is much smaller 
a p(l/ i l) = 2^~g2^2] — p a e(|Ml) = ~2~B^R^\ — f ^ lan tne thickness d of the blast wave region we will establish 
p ^ ' that the inflowing proton-electron-beam is effectively isotro- 
Integrating Eqs. ([33]) and (|4|) over t and fi gives because of pised in the blast wave plasma. 
a p , e (£t = -1) = Inserting Eqs. @ and @ we find with 

* , , , r-l/Rpk i i 

dfi [F PiC (fi ,t)-F p , ( M ,t = 0)] / dfX [ Fp ( fX> T i )-F p {n,t = 0)} = --(l + —)(42) 



-l J-i z -n-P 

and 

1/R c k 



*p^(H)5Z / fiJiGiW*!- 1 ,*)-^ (36) 



respectively. Evaluating Eq. ( ^6| ) at /i = —(R p k) 1 for protons 
and at /i = (i? c fc) _1 we obtain 



y 3 W(^,77)-^c(M = 0)] = -i(l--^) (43) 



that Eq. ( [39[ ) becomes 

/•t ap 



0^ 1 (|i2 p fc|- 1 ) / ' ^V^pGM) - -Fp(M = 0)] (37) x £ + ^[Ifcl-^x 1 ] 

7-1 x=o,p 
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m c Va ^ p ,c B 2 R ™i 
m„ 2v c 2 k 2 p nbT 



x J2 [l + ikR^mikl-R- 1 ] 

x— e,p 

~ -Z k- 2 (44) 
The general solution of Eqs. (^0|) and (^T]) at time t = T± is 

= yjY + ${Z + I + {0)-I-(0))* 

+0.5 (Z + I+(0) -J_(0)) (45) 

= ^ + i(^ + /+(0)-/-(0))2 

-0.5 (Z + 7+(0)-Z_(0)) (46) 

where 

F = I+(0)7_(0) (47) 

If the initial turbulence is much weaker than the self-generated 
turbulence I{k, 0) << |Z(fc)| and has a vanishing cross-helici- 
ty/+(fc,0) =I-(k,0) =J(fc,0) we obtain for Eqs. (g5j||0) 
approximately 



(48) 



According to Eq. (^4|) Z(fc) is negative so that Eq. (^8j) reduces 
to 

l 2 (k,0) 



/+(fc ' Tl) " |Z(fc)| 
and 

I-(k,Ti)~ \Z(k)\ 



(49) 



(50) 



i.e. the beam mainly generates backward moving Alfven waves 
in the blast wave plasma. 

2.3. Energy budget 

The total enhancement in magnetic field fluctuation power due 
to proton and electron isotropisation is obtained by integrating 
Eq. @ using Eq. (Q 



dk[l-(k,Ti) - l-(k,0)\ 

X) 

dkZ{\k\) 

— OO 

m p vc 2 «bT m p 

47rV r An i mpur[l + (— )] 



(51) 



so that the change in the magnetic field fluctuation energy den 
sity is 

(SB-) 2 1 



8tt 2 



V Ani m p vT[l + (-^)] 



(52) 



Alfven waves possess equipartition of wave energy density be- 
tween magnetic and plasma velocity fluctuations, so that the 
total change in fluctuation energy density due to pitch angle 
isotropisation is 



All = AU Sv + AU SB = 2AU SB = V A nim p vT[l - 



(53) 



For consistency we show that this increase in the energy density 
of the fluctuations is balanced by a corresponding decrease in 
the energy density of the beam protons and electrons during 
thei r isotr opization. We follow here the argument of Bogdan et 
al. (1991), made originally for pick-up ions in the solar wind, 
and generalise it to relativistic beam velocities. 

As the beam particles scatter away from their initial pitch 
angle a — ir and speed V, at each intermediate /j, = pii jp they 
are confined to scatter approximately on a sphere centered on 
the average wave speed of waves which resonate with that pitch 
angle, given by (Skilling 1975, Schlickeiser 1989) 



V ph = V A 



i + {(r^)- 1 ) - i-dR^r 1 ) 

/ + ((i2 i M)- 1 )+7-((^M)- 1 ) 



-V A 



(54) 



because according to Eqs. (^9j|50|) the intensity of backward 
moving waves is much larger than the intensity of forward 
moving waves. Thus the beam particles scatter in momen- 
tum space on average onto a surface axisymmetric about e z 
bounded by a sphere centered at — Va, and passing through the 
ring of beam particle injection. Once the beam particles are uni- 
formly distributed over the surface Vf(fi), their energy density 
wt is 



d 3 pf t {p, A*, t = U) -fiiriiC 2 



(55) 



where the surface Vf(fi) is the relativistic sum of the initial 
velocity V and 



V 2 = I d f xV ph ( l i)=-V A ( l i + l) 

' cos(a) 



We readily find 

V + V 2 



(56) 



(57) 



or in units of the speed of light c (f3 = V/c, Pf(p) — Vf/c, 
Pa — Va/c) we obtain 



AM 



(58) 



1 - pp A ([x + 1) 

The phase space distribution function of the beam particles af- 
ter isotropisation is 



fi(p,fx,t = ti) = S (p-Pt(p)) 

with pi — miC/3fjf and 7f = (1 — p 2 )^ 1 / 2 . Inserting Eq. 
in Eq. (p5|) we obtain 



(59) 



Wi 



E 



n\raiC 



2 /■! 



d^di) 



Martin Pohl and Reinhard Schlickeiser: Channelled blast wave model for AGN 



7 



E 



•n\m,iC 



dp 



i - pp A (n + 1) 



^[i-pp A {^ + \) 2 -[p-P(^ + \)] 2 

Substituting u = p + 1 we derive 

- 2 r2 i pp a 



(60) 



E 

i 

E 



du- 



y/l - P\u 2 



arcsin(2/? A ) - P + P^fl 



(61) 



For our magnetic field strength B ~ 1 G and blast wave densi- 
ties ?ib — 10 6 n b 8 /rioo cm -3 , the Alfven speed p A — 10~ 2 is 
much less than unity, so that we may expand Eq. ( j6l| ) for small 
values of p A yielding 



Wl 



n imi c 2 r [l + -pl + -p% - pp A (i - pi) 

i 

J2n imi c 2 r[i-p A p] 



(62) 



The initial energy density of the beam particles wq can be cal 
culated using the initial beam distribution function (|l|), i.e. 

n 



fi(p, fi,t = Q) = ^S(p + l)S(p - P) 
yielding 

w = E / ^ ^' 1 = Q)"fimiC 2 = E n i m »c 2r (64) 

Obviously the change in the energy density of the beam parti- 
cles then is 

Aw = Wl — Wq = — ni1TliC 2 TP A P 

i 

TYl 

= - ni m p vV A T[l + — ] (65) 



which is exactly —AU from Eq. (p3[). The plasma turbulence 
is generated at the expense of the beam particles which relax 
to a state of lower energy density; or with other words, the ex- 
cess energy density in the non-isotropic beam distribution is 
transferred to magnetohydrodynamic plasma waves that scat- 
ter the beam particles to an isotropic distribution in the blast 
wave plasma. If this isotropisation is quick enough - which 
we will calculate next - the beam particles attain an almost 
perfect isotropic distribution function ( |59| ) in the blast wave 
plasma, because to lowest order in p A the final particle mo- 
mentum pf ~ PTrriiC is independent of p. 

2.4. Isotropisation length 

Using the fully-developed turbulence spectra ( (49[pT)| ) in Eq. 
(18) we obtain for the pitch angle Fokker-Planck coefficient 
of the beam particle i 



D„ 



nn 2 1 



/' 



B 2 



(66) 



where 



Z ^ 7 i 



r * V VII — 



1 Va ) 

\hvtV a \ 



I- 1 



vp,+V A 



1 \iiv-V A \ 
\Z \\vp + V A \ 



\pv + V A \ 

\Z \\vp-V A \ 3 



(67) 



For ease of exposition we assume that the initial turbulence 
spectrum has the form I(k, 0) = Iok~ 2 , so that Eq. ( |67| ) be- 
comes 



\Zo\l 

n 2 



\vp + V A \+ A\vp-V A \ 



(68) 



where A = (I /\Z \) 2 . 

According to quasilinear theory (e.g. Earl 1973, Schlick- 
eiser 1989) as a consequence of pitch-angle scattering the beam 



particles adjusts to the isotropic distribution < J4 1| ) on a length 
scale given by the scattering length 



A = ^ 
8 



-1 MM 



3 B 2 



(/i) 47T |Z„| 



J(e,A) 



(63) with e = V A /v and 



J(e,A) = / ^(l- M 2 )[| M + e| + A|/x-ei 



(69) 



(70) 



To obtain Eq. (|69|) we have inserted Eqs. ( |66| ) and (|58j). After 
staightforward but tedious integration the value of the integral 
( p70| ) to lowest order in the small parameters e < < 1 and A < < 
1 is 

J(e, A) ~ -21n(v / 2eA) (71) 
so that the scattering length (59) becomes with Eq. (Q) 

3 B 2 



A 



2^ |Z | 
3 m D 



In(v^cA)] 
2 w nbF 



7T TO e Wp iC i?p Va ™i 

= l^^jn^eA)" 1 



■ln(\/2eA)] 



Inserting our typical parameter values we obtain 

1/2 

A = 2.2 • 10 11 — \n{\/2p A A)~ 1 cm 



(72) 



(73) 



Note that the ratio of initial to fully developed turbulence in- 
tensities enters only weakly via the logarithm. For turbulence 
intensity ratios from 10 _1 to 10~ 5 implying values of A from 
10~ 2 to 10~ 10 we find with A ~ 10~ 2 that ln[V2p A A]^ 
varies between 8.9 and 27.3. Taking the larger value in Eq. ( J73[ ) 
yields for the scattering length in the blast wave plasma 



A ~ 6- 10 



in 



1/2 
n b,S 



(74) 
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which corresponds to an isotropisation time scale of 



t R = X/c 



i 1/2 

'b,8 



(75) 



If the thickness d of the blast wave region is larger than the 
scattering length (75), indeed an isotropic distribution of the 
inflowing interstellar protons and electrons with Lorentz factor 
(see Eq. @) < Y >= Y{1 - (3 A (3) ~ T in the blast wave 
frame is effectively generated. In the following sections we in- 
vestigate the radiation products resulting from the inelastic in- 
teractions of these primary particles with the cold blast wave 
plasma. 

This discussion of particle isotropization on self-excited 
turbulence applies to both AGN and GRBs. The radiation mod- 
elling of GRBs often requi res en ergy equipartition between 
electrons and protons (Katz 1994), based on detailed plasma 



physics considerations (e.g. Beloborodov & Demiahski 1995 



Smolsky & Usov |1996fc Smolsky & Usov |1999| ), However, the 
isotropization itself provides electrons with an energy roughly 
2000 times smaller than that of the protons. Therefore a reac- 
celeration of electrons would be required to reach equipartition. 
Since the energy density of the turbulence is only a fraction 
(oc (3a) of the energy density of the incoming beam, the tur- 
bulence is not energetic enough to reaccelerate electrons for 
Pa <C 1, and further studies may be required to understand 
the production of radiation in GRBs. The situation is differ- 
ent with AGN, for which equipartition is not required and for 
which the variability timescales are substantially longer. Here 
we will discuss the radiation modelling for AGN on the ba- 
sis of the particle distributions resulting from the isotropization 
process alone. 

3. Radiation modelling of the blast wave 

Now we deal with isotropic particle distribution functions in 
the blast wave frame, which itself is not stationary because the 
blast wave sweeps up matter and thus momentum. Momentum 
conservation then requires a deceleration of the blast wave de- 
pending on whether or not the swept-up particles maintain their 
initial kinetic energy, and depending on possible momentum 
loss from anisotropic emission of electromagnetic radiation. 

As we have sketched in Fig.[l], the blast wave is assumed to 
have a disk-like geometry with constant radius R and thickness 
d that moves with bulk Lorentz factor T. The matter density in 
that disk n b is supposed to be orders of magnitude higher than 
that of the ambient medium n* . 

In the blast wave frame the external density n, = Y n* and 
sweep-up occurs at a rate 



iV( 7 ) = nR 2 cn* \/Y 2 - 1 5{j-Y) . 



(76) 



As we have seen in the preceding section, the relativistic elec- 
trons and protons get very quickly isotropised and lose little 
energy in the process. This has two consequences: the sweep- 
up is a source of isotropic, quasi-monoenergetic protons and 
electrons with Lorentz factor Y in the blast wave frame. The 



isotropisation also provides a momentum transfer from the am- 
bient medium to the blast wave. 



3.1. The equation of motion of the blast wave 

In a time interval St the blast wave sweeps up a momentum of 

SU = 7T R 2 rnp c 2 n* (r 2 - 1) St (77) 

which is transferred from the swept-up particles to the whole 
system with mass 

M B w = ttR 2 d (n b rn p + 7i nth m nt h) (78) 
where the density and relativistic mass of the energetic particles 

(79) 



«nth TOnth = "lp / d7 711(7) 



have to be added to that of the thermal plasma. Therefore, the 
blast wave will tend to move backwards and its Lorentz factor 
relative to the ambient medium is reduced to 



r' = rwi 



/ 511 



\ M B w c 



- Vr 2 -! 



sn 



M B w c 



(80) 



which can easily be integrated numerically. For a highly rela- 
tivistic blast wave we can expand the blast wave equation of 
motion and derive the timescale for slowing down 



6V_ 

St 



Yd (n b m p + n nth m nth ) 
n* m p c (r 2 — l) 3 / 2 



(81) 



In the laboratory frame this timescale is longer by a factor Y. 

We may also integrate the inverse of the slow-down rate to 
derive the observed time frame of the deceleration, that is the 
time at which the blast wave would be observed with a particu- 
lar Lorentz factor. This may be interesting to compare with the 
results of VLBI observations of blazars. In the radiative regime, 
i.e. when the internal energy is radiated away quickly, we may 
neglect the mass loading for Y ^> 1 and obtain 



rib d 
nf c 



Si 



dY 



, r'(i-/3V) 
(r 2 - 1) 3 / 2 



(82) 



where we have assumed an interstellar medium with a constant 
density n*. In Fig. ^ we show typical solutions of the inte- 
gral Eq.(p2[). The observed time needed for a deceleration to 
Lorentz factors Y of a few is independent of the initial Lorentz 
factor Yq 3> T, and it varies much less strongly with the aspect 
angle than did the initial Doppler factor. 

We can also calculate the distance traveled by the blast 
wave. For this we replace the retardation factor in Eq.(|82|) by 
(3c and we obtain for Yq 3> Y 



L 



n h d Y 
In — 

2n* T 



(83) 
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10 100 
Lorentz factor T 

Fig. 2. The observed time plotted versus the Lorentz factor of 
the blast wave. The solid line is calculated for an observer lo- 
cated in the direction of the jet, i.e. (i=l t and initial Lorentz 
factor 300. The dotted line, derived for an initial Lorentz fac- 
tor of 100, and the dashed line, calculated for an aspect cosine 
fi=0.9, show that the observed time needed for the decelera- 
tion to semirelativistic velocities is independent of the initial 
Lorentz factor and varies little with the aspect angle. 



Comparing the swept-up mass with the initial mass of the blast 
wave we see that mass loading in the radiative regime is im- 
portant only for V & 2, and thus its neglect in the derivation of 
Eqs. |82| and [33] is justified. If the blast wave is strong enough, 
L will be in the range of a few kpc and the blast wave may 
leave the host galaxy at relativistic speed. It may travel through 
the dilute intergalactic medium for several Mpc, before it even- 
tually runs into a denser gas structure and converts its kinetic 
energy into radiation. Thus the power of the blast wave, i.e. the 
luminosity, is linked to the probability of escape from the host 
galaxy. When the systems are viewed from the side, they would 
fit into the morphological classification of Fanaroff-Riley I and 
Fanaroff-Riley II galaxies: while the former are less luminous, 
they are also core dominated, whereas the latter exhibit very 
luminous radio lobes. 



3.2. The evolution of the particle spectra 

Eq. ( f76| ) states the differential injection of relativistic particles 
in the blast wave. Here we concentrate on the protons because 
they receive a factor of m p /m e more power than electrons, 
which also have a low radiation efficiency for 7 <C 1000. Elec- 
trons (and positrons) are supplied much more efficiently as sec- 
ondary particles following inelastic collisions of the relativis- 



tic protons. Since no reacceleration is assumed, the continuity 
equations for protons and secondary electrons read 



h — — 



dt 



7^(7) 



JV P (7) , ^p(7) 



T 



JVp( 7 ) (84) 



N e ("f) (85) 



where Te is the timescale for diffusive escape, and Tn is the 
loss timescale associated with (p — ► n) reactions with neutron 
decay occurring outside the blast wave. For positrons catas- 
trophic annihilation losses on a timescale T ann are taken into 
account. The injection rate of secondary electron is related to 
the rate of inelastic collisions. 

Energy losses arise from elastic scattering at a level of 



7el,p 



3-10" 



-16 



rib 



7 



Vl 2 



1 



(86) 



where we have assumed the ambient medium as a mixture of 
90 % hydrogen and 10 % helium. The energy losses caused by 
inelastic collisions can be determined by integration over the 
secondary particle yield. Here we use the Monte Carlo model 
DTUNUC (V2.2) (Mohring & Ranft 
Ferrari et al 



1996a 



1991 



1997) 



1994 



Ranft et al 
which is based on a 
1994b. This MC model for 



Engel et al. 
dual parton model (Capella et al 
hadron-nucleus and nucleus-nucleus interactions includes var- 
ious modern aspects of high-energy physics and has been suc- 
cessfully applied to the description of had ron production in 
high- energy collision s (Fer rari et al. 1996b , Ranft & Roesler 



1994, Mohring et al. 1993 Roesler et al. 1998). The total en 



ergy losses from inelastic collisions are well approximated by 



' 7inel = 7 



7 + 1 



The timescale for neutron escape after (p 
relativistic energies is approximately 



3-10 



15 nr 1 



exp 



L 



7 2.7- 10 13 cm 



(87) 



n) reactions at 



sec (88) 



where L is the neutron path length through the blast wave and 
the brackets denote an average over the neutron emission angle. 
In all the examples discussed in this paper the exponential in 
Eq. ( |8~8| ) will be close to unity. 

The Monte-Carlo code also provides the differential cross 
sections for the pion production, on which we base our calcu- 
lation of the pion source functions. Neutral pions decay imme- 
diately into two 7-rays and the 7-ray source function 



/•OO 

97(f) - 2 / 



El- 



(89) 



Charged pions decay into muons and finally into electrons or 
positrons. The secondary electron source functions are essen- 
tially determined by the kinematic of the two decay processes 
(Jones 1963 Pohl 1994). The secondary electrons loose energy 
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mainly by inverse Compton scattering, synchrotron emission, 
bremsstrahlung, and elastic scattering. 



C/ph 



7ic = 2.7- HT 14 

-7s y = 2.7-10 
- 7Br = 8 • 1CT 16 n h ( 7 - 7- 1 ) 



14 u b 



(7 2 



"7el,i 



6-10 



-13 



7 



(90) 

(91) 

(92) 
(93) 



where the energy densities Ub and U p h are in units of eV/cm 3 . 
The pro bability of annihilation per time interval is (Jauch & 
Rohrlich 11976k 



T: 



3 c er T rib 



!(7+l)v / 7 T ^l 

( 7 + 4 + 7 -i) hi (7 + V7 5 



(94) 



/8(7 + 3) 



where we have assumed that the temperature of the background 
plasmas is higher than about 100 eV, so that a Coulomb cor- 
rection and contributions from radiative recombination can be 
neglected. 

The diffusive escape of particles from the blast wave scales 
with the scattering length A as derived in Eq.(f74|). For the es- 
cape timescale we can write 



T ¥ . = — = 



3d 2 



D \(3c 



1.67- 10~ 19 



j3 Jnh 



sec 



(95) 



If escape is more efficient than the energy losses via pion pro- 
duction, then the bolometric luminosity of the blast wave will 
be reduced by a factor Te/t w ~ Te |7mel|/7- 

The proton spectra evolve differently depending on the rel- 
ative effect of the losses and the blast wave slow down. If there 
were no losses, the relativistic mass loading Eq. ( |79[ ) would be 
trivial to calculate. Combining the mass loading equation with 
the deceleration equation Eq. ( |80[ ) we obtain after some calcu- 
lus in the limit T > 1 



r = 



c r 4 

n h dT 



(96) 



where Tq is the initial Lorentz factor of the blast wave. By com- 
parison with the simple timescale Eq. (^Tj) we see that the mass 
loading reduces the deceleration by a factor T/Tq. The proton 
continuity equation (Q) can be trivially integrated in the ab- 
sence of losses to 



AT( 7 ) 



drNpij) =C 



dT 



p/-4 



c'r 



for r < 7 < r 



(97) 



So if the energy losses can be neglected, the blast wave deceler- 
ation alone would cause approximately a 7" 3 spectrum in the 
relativistic range. 



3.3. Optical depth at - -ray energies 

Here we discuss a system of rather peculiar geometry, a very 
thin disk, and a strong angular variation of optical depth ef- 
fects can be expected. While for photons traveling exactly in 
the blast wave plane the line-of-sight through the system has an 
average length R, and the system may appear optically thick, 
the angle-averaged line-of sight has a length 



Z=^ln 
2 



2R 
~~d 



(98) 



which is approximately equal to d. Therefore the total optical 
depth may be small because the photons would interact only in 
a very small solid angle element. In all the cases presented here 
this is true for the Thompson optical depth t = a^n c L, so that 
Compton scattering will not effectively influence the emission, 
though <tt n c R may be of the order of 1 . 

At photon energies above 511 ke V in the blast wave frame 
photon-photon pair production may occur and inhibit the es- 
cape of 7-rays from the system. This process is not important 
in our model, because very high apparent luminosities can be 
produced with moderate photon densities in the source as a re- 
sult of the highly relativistic motion of the blast wave and the 
associated high Doppler factors. For example, given a system 
radius R = 10 15 cm and a Doppler factor D = 100, an appar- 
ent luminosity L = 10 48 erg sec -1 at MeV energies in the blast 
wave rest frame corresponds to an optical depth of r 77 ~ 10~ 4 . 



4. High energy emission from the blast wave 

In Fig. U we show the reaction and radiation channels of the 
swept-up particles together with the approximate photon en- 
ergies for the emission processes. Since the source power of 
secondary electrons and the power emitted in the form of 7r°- 
decay 7-rays is the same, the bulk of the bolometric luminosity 
will be emitted in the GeV to TeV energy range, independent 
of the choice of parameters. Thus the unexpected finding, that 
7-ray emission of blazars generally provides a major part of the 
bolometric luminosity, would find a natural explanation if as in 
this model pion production in inelastic collisions were the main 
source of energetic electrons. 

In this paper we concentrate on the high-energy 7-ray emis- 
sion, and especially on the processes of ir° decay, bremsstrah- 
lung, and pair annihilation, for they scale only with internal 
parameters like the gas density rib- A discussion of inverse 
Compton scattering and neutrino emission will be presented 
elsewhere. The inverse Compton process must inevitably oc- 
cur, but its rate, the spectral and angular distribution depend on 
the actual choice of target photon field. If this scattering pro- 
vides significant energy losses for electrons and positrons, the 
solution of Eq. (|5]) will be particularly arduous. 

The 7T° production spectrum is calculated using the Monte- 
Carlo code DTUNUC described earlier in this paper. 
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E p= r o m P c 2 
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Fig. 3. Scheme of the reaction and radiation 
channels of relativistic blast waves. Here the 
logical flow is from top to bottom. In each 
box we also note the approximate energy of 
the particles and the photons in the blast 
wave frame. Depending on the Doppler factor, 
i.e. on the aspect angle, the observed energy 
may be substantially higher. The channels dis- 
played in solid line boxes will be discussed in 
detail in this paper, whereas those printed in 
dashed line boxes will be the subject of forth- 
coming publications. 



The bremsstrahlung spectrum is determined using the rel- 
ativistic limit of the differential cross section (Blumenthal & 



Gould 1970) 



da 

de bs 
where 



e 



!> u = 4 In 



4 e 
37 

27 



(99) 

- 2 

and the photon energy e is in units of m e c 2 . The annihilation 
spectrum is obtained from the differential cross section (Jauch 
& Rohrlich 11976b 



1-1 



da 
~dl 

with 



AN 



P 2 ! 2 



e 
ci 



e 



2(1 + 7) A + 7 



(100) 



e + ej = 1 + 7 and 



^ T ^r<-<7(l + /3) 
7(1 + 0) ei 



We follow the evolution of the electron and positron spectra to 
Lorentz factors of 3, below which the bremsstrahlung emissiv- 
ity is not calculated and the annihilation spectrum is derived in 
the non-relativistic limit. 

Synchrotron emission can be expected in the optical to X- 
ray frequency range. For head-on jets synchotron emission may 
be observable at X-ray energies. The peak of the observed syn- 
chrotron spectra in vF v representation scales roughly as 



V l lx " ( Gauss) (200) (2OO 



kcV 



(101) 



with D denoting the Doppler factor, but free-free absorption 
and the Razin-Tsytovich effect, in addition to synchrotron self- 
absorption, would inhibit strong synchrotron emission in the 
near infrared and at lower frequencies. Radio emission may 
become observable at later phases, when the blast wave has de- 
celerated and the Doppler factor is reduced, or when the blast 
wave medium is diluted, e.g. as a result of imperfect collima- 
tion. 

The peak of the observed 7r°-decay spectrum in vF v repre- 
sentation scales as 



f- 

V 100 



D 
100 



TeV 



(102) 



Comparing Eqs.( J101| 9 and ( |102[ ) we see that hard X-ray syn- 
chrotron emission as appare ntly o bserved from the BL Lacertae 
object Mkn 501 (Pian et al. |1998| ) implies that the vF v peak of 
the high energy emission is located far in the TeV range of the 
spectrum. At this point we like to issue a warning to the reader, 
that one has to be careful when comparing the model spectra 
to actual data. The TeV 7-ray spectra of real sources, even the 
closeby ones, are modulated during the passage through the 
intergalactic medium, for the 7-ray photons undergo pair pro- 
duction in collision with infrared background photons. Unfor- 
tunately the intergalactic infrared background spectra, and thus 
the magnitude and the energy dependence of the optical depth, 
are not well known. The current limits are compatible with se- 
vere absorption at all energies above 1 TeV even for Mkn 421 
and Mkn 501. Therefore, model spectra which display a vF v 
peak at 10 TeV are not incompatible with the observed peak 
energies in the range of ~ 0.5 TeV for Mkn 421 and Mkn 501. 
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Fig. 4. Spectral evolution of a relativistic blast wave in an environment of constant density. The top row displays F v spectra of the 
7-ray emission, whereas the bottom row shows vF v spectra from the optical to high energy 7-rays including also the synchrotron 
emission. The solid lines display 7r°-decay 7-rays, the dot-dashed lines bremsstrahlung, the dashed lines annihilation emission, 
and the dotted lines represent the synchrotron emission. From left to right the panels show the spectra after one hour, 10 hours, 
and 100 hours observed time. The parameters are: ro=300, d=3 • 10 13 cm, R=10 14 cm, B= 2.0 G, n*=0.2 cm~ 3 , n^=5 ■ 10 8 cm -3 , 
for an AGN at z=0.5 viewed at an angle (9 b s = 0.1°. After 100 hours the blast wave has decelerated to T ~ 106. 



In Fig. ^ we show the spectral evolution of high energy 
emission from a collimated blast wave for a homogeneous ex- 
ternal medium. Even for very moderate ambient gas densities 
the high energy emission will be very intense. Some general 
characteristics of our model are visible in this figure. The tt°- 
decay component dominates the bolometric luminosity. This is 
a direct consequence of the hadronic origin of emission, as the 
source power available to the leptonic emission processes is al- 
ways less than the pion source power, for the neutrinos carry 
away part of the energy. 

The energy loss time scale of the high energy electrons is 
always smaller than that of the protons, and therefore the X- 
ray synchrotron intensity will follow dlosely variations in the 
TeV 7r°-decay emission. Thus there should be a general rela- 
tion between the X-ray and the TeV appearance of sources with 
jets pointing towards us. Note, however, that an observer may 
see an imperfect correlation if a) the synchrotron and the ir°- 
decay spectra are not observed near their i/i^-peaks, or b) in- 
verse Compton scattering contributes significantly to the X-ray 
emission, or c) strong absorption shifts the apparent ^F^-peak 
of the 7r°-decay emission. Correlations between the X-ray and 
TeV 7-ray emission of BL Lacertae objects have indeed been 
observed in case of Mkn 421 (Buckley et al. 1996 1 and Mkn 
501 (Catanese et al. 1997, Aharonian et al. 1999 and Djannati- 
Atal et al. 11999b, but at least for Mkn 421 recent data show that 



there is no one-to-one correlation during flares (Catanese et al. 



1999). 



The high Lorentz factors in our model imply that the ob- 
served emission depends more strongly on the aspect angle 
than in conventional jet models with V ~10. This is shown in 
Fig. ||, which displays the spectral evolution of a source viewed 
at an aspect angle of 2° (/i=0. 99939). The model parameters 
have been changed slightly compared with those used for the 
head-on case (Fig. Q), mostly to make the source brighter. 
The fundamental difference to the 9 ~ 0° case is, however, 
that the Doppler factor D increases with decreasing T when 
/1 > f3. In contrast to Fig. f|, where always D > V, here 
D increases from an initial value of ~ 5 to a maximum of 
D max = (1 — /i 2 )~ 5 ~ 29. As a consequence the apparent 
evolution of the source spectrum is much slower, and the peak 
energies are smaller and don't change much during the evolu- 
tion. The synchrotron peaks are located in the optical, and the 
x-ray luminosity is small. 

Strong TeV emission implies that the systems need to be 
observed head-on. Even for T = 200 an aspect angle 9 < 1° is 
required. Viewed at aspect angles of a few degrees the sources 
look like the typical hard spectrum EGRET sources: the spec- 
trum peaks somewhere in the GeV region, the sources reach 
their maximum luminosity when the blast wave has decelerated 
to Lorentz factors around ~ sin -1 9, the luminosity spectrum 
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Fig. 5. The same system as in Fig. ^ viewed by an observer located — 2° off the blast wave direction of motion (/i=0. 99939). 
Note that now the panels show the spectra after 10 hours, 100 hours, and 1000 hours of observed time. Here also d=10 14 cm, 
R=2 • 10 15 cm, n*=1.5 cm~ 3 , and rib=10 8 cm -3 to make the source appear as bright as before. The system appears to evolve 
much slower than for 8 — 0.1°. Note that the left column here represents the same observed time as the middle column in Fig. || 
The spectral evolution differs also qualitatively from that for 9 — 0.1°, for the blast wave deceleration causes an increase of the 
Doppler factor as long as /i < j3, which in this example holds as long as the observed time t< 7 • 10 6 sec. 



displays a minimum in the keV to MeV range, and the temporal 
evolution is slower. 

The dependence of the vF v peak of the high energy 7-ray 



emission on the aspect angle Eq. (102) also implies a general 



relation between the observed E„o max and the a ppare nt super- 



luminal velocities in VLBI images. Using Eq. ( |102[ ) we may 
write 
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From this limit we derive for the apparent velocity 
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This is in accord with the observation that detected EGRET 
sources, i.e. primarily MeV - GeV peaked emitters, often show 
rapid superluminal motion, whereas the only TeV peaked 7-ray 
source with good VLBI coverage, namely Mkn 421, displayed 
subluminal motion in the time range of 1994 to 1997 (Piner et 
al. p999b. 



All spectra exhibit a strong annihilation line at an energy 
of D ■ 511keV. The copiously produced positrons cool down 
to thermal energies before annihilating, and a relatively narrow 
line is produced at a given time. However, since the Doppler 
factor may vary during the integration time of current 7-ray 
detectors (typically three weeks (!) in case of COMPTEL), the 
line may look like a broad hump in the data. 

The spectra shown so far illustrate the effect of the blast 
wave deceleration and the thus modulated injection rate. To 
this we oppose in Fig. || and Fig. (7] the case of the blast wave 
interacting with an isolated gas cloud. The light curves of high 
energy 7-ray sources often display phases of activity lasting for 
months or years, on wh ich rapid fluctuations are superimposed 
(e.g. Quinn et al. 1999 ). We propose that the secular variabil- 
ity is related to the existence or non-existence of a relativistic 
blast wave in the sources, and thus to the availability of free en- 
ergy in the system. Strong emission on the other hand requires 
that the available kinetic energy of the blast wave is converted 
into radiation. Depending on the distribution of the interstel- 
lar matter the kinetic energy of the blast wave may be tapped 
only every now and then. The observed fast variability of high- 
energy 7-ray sources would thus be caused by density inhomo- 
geneities in the interstellar medium of the sources. The situa- 
tion displayed in Fig. ^ would then correspond to one of the 
rapid outbursts in the 7-ray light curves of AGN, which would 
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LOG(E [GeV]) LOG(E [GeV]) LOG(E [GeV]) 

Fig. 6. Spectral evolution of a relativistic blast wave having traversed a gas cloud of density n*=0.2 cm~ 3 and thickness 5 • 10 16 
cm. The other parameters are: ro=300, d=10 14 cm, R=3 • 10 14 cm, B= 2.0 G, nb=10 9 cm~ 3 , for an AGN at z=0.5 viewed at 
an angle — 0.1°. The panels show from left to right the spectra after 0.05 hours, 0.5 hours, and 5 hours observed time. The 
Lorentz factor of the blast wave did virtually not change, hence repeated cloud crossings would produce multiple outbursts with 
the same spectral evolution. 




LOG(E [GeV]) LOG(E [GeV]) LOG(E [GeV]) 

Fig. 7. The same system as in Fig. || viewed by an observer located 2°off the blast wave direction of motion (/.t=0. 99939). Note 
that now the panels show the spectra after 0.5 hours, 5 hours, and 50 hours of observed time. Here also d=3 • 10 14 cm, R=4 • 10 15 
cm, n*=l cm -3 , and T=30 to make the source appear as bright as before. The choice of the initial Lorentz factor was done to 
maximize the observable flux. 
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Fig. 8. Evolution of the TeV spectra for a single outburst with 
similar parameters as in Fig.|| The curves represent the spectra 
at different times after the beginning of the cloud crossing. The 
labels refer to observed time. 
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generally follow the gas density variations in the volume tra- 
versed by the blast wave. 

The large collection area of currently operating Cerenkov 
telescopes permits studies of spectral variability on very short 
timescales. The observational situation is disturbing, however, 
since data taken with different telescopes can be seemingly in- 
consistent. As an example, both the CAT team and the HEGRA 
group performed observation of Mkn 501 in the time of March 
1997 to October 1997. The CAT data show a statistically sig- 
nificant correlation between the spectral hardness and flux in 
the energy range between 0.33 TeV and 13 TeV. The HEGRA 
data on the other hand give no evidence in favor of a correlation 
in the energy range between 1 TeV and 10 TeV. 

Both behaviours could be reproduced with our model. In 
Fig. H we show the spectral variability at TeV energies only. 
This plot is essentially a blow-up of Fig. ^. It is obvious that the 
rise time of an outburst is much smaller than the fall time. Also 
in the declining phase the source undergoes strong spectral evo- 
lution, a consequence of the proton energy loss timescale being 
smaller than the escape timescale in this example. Would the 
source repeat the outbursts every, say, five hours, we would def- 
initely observe hard flare spectra and softer quiet phase spectra. 

Were the clouds bigger, so that the cloud crossing time 
in the blast wave frame is longer than the proton energy loss 
timescale, and would we allow for a dilute intercloud medium, 
both during the passage through the cloud and while travers- 
ing the intercloud medium the TeV spectra would approach 
a steady-state. If we further chose the parameters such that 
the blast wave deceleration is small, we would reproduce the 



Fig. 9. Evolution of the TeV spectra for multiple outbursts re- 
sulting from clouds embedded in the interstellar medium. The 
interstellar density n* was assumed as 10~ 2 cm -3 in the inter- 
cloud medium and 0.5 cm~ 3 in the clouds, which were sup- 
posed to be 1 pc in size with a cloud-to-cloud distance of 
10 pc. Here also T=100, d=10 14 cm, R=3 • 10 14 cm, n h = 
3 • 10 9 cm- 3 , 6 = 0.1°, B=2.0 G, and z=0.5. The modeled 
data have been divided into two groups according to the flux 
above 0.3 TeV: the flare spectra, here shown as solid line, and 
the quiet spectra, plotted as a dotted line. The top panel show 
the ratio of the two, which is compatible with a constant for en- 
ergies above 0.3 TeV, thus reproducing the HEGRA result for 
Mkn 501 (Aharonian et al. 



19991 



HEGRA result of variability without spectral evolution in the 
TeV emission from Mkn 501. This is shown in Fig. ^, where 
the modeled data have been divided into two groups according 
to the flux above 0.3 TeV: the flare spectra and the quiet spec- 
tra. This is the same classification that was used in the HEGRA 
data analysis, except their use of five groups instead of two. It 
is obvious that the spectra in both groups are virtually identical 
in the energy range that was used for the classification. 

5. Conclusions 

We have shown that a relativistic blast wave can sweep-up am- 
bient matter via a two-stream instability which provides rela- 
tivistic particles in the blast wave without requiring any accel- 
eration process. We further have demonstrated that the blast 
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waves can have a deceleration timescale that allows escape 
from the host galaxy, and hence allows the formation of gi- 
ant radio lobes without invoking additional energisation. While 
the relativistic blast wave is rather long-lived, the swept-up rel- 
ativistic particles in it have short radiation lifetimes. The abun- 
dance of such particle, and hence the intensity of 7-ray emis- 
sion, traces the density profile of the traversed interstellar gas. 
We imagine that the secular variability of the 7-ray emission 
of AGN is related to the existence or non-existence of a rela- 
tivistic blast wave in the sources, and thus to the availability of 
free energy in the system. The observed fast variability on the 
other hand would be caused by density inhomogeneities in the 
interstellar medium of the sources. Since we do not consider 
any re-acceleration of particles in the blast wave, the evolution 
of particles and the blast wave is completely determined by the 
initial conditions. 

The high energy emission of a relativistic blast wave mov- 
ing approximately towards the observer has characteristics typ- 
ical of BL Lacertae objects. In particular, 
- the high energy spectra are very hard with photon indices 
< 2, in accord with the unspectacular appearance of TeV-bright 



sources at GeV energies (Buckley et al. 1996) 
- as can be seen in Fig. ^, observable increase and decrease 
of intensity at TeV energies can be produced on sub-hour time 
scales, in accord with the observed variability time scales of 
Mkn 421 (Gaidos et al. 



1996) 



- for multiple outbursts the intensity can follow the variation of 
the ambient gas density with little spectral variation, simila r to 
the observed behaviour of Mkn 501 (Aharonian et al. 



1999) 



x-ray synchrotron emission is produced in parallel to the 7- 



rays as was observed from Mkn 501 (Pian et al. 1998) 
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